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^ . We prove that for the Ising model on a lattice of dimensionality d > 2, the zeros of the 

■ partition function Z in the complex /x plane (where /i = e~^^^) lie on the unit circle = 1 

^ ■ for a wider range of Knn' = /3<^nn' than the range Knn' > assumed in the premise of the 
J-^ ■ Yang-Lee circle theorem. This range includes complex temperatures, and we show that it 
I is lattice-dependent. Our results thus complement the Yang-Lee theorem, which applies for 
^ I any d and any lattice if J„„' > 0. For the case of uniform couplings Knn' = K, we show 
', that these zeros lie on the unit circle = 1 not just for the Yang-Lee range < m < 1, but 
also for (i) —Uc,sg < m < on the square lattice, and (ii) —Uc,t < u < on the triangular 
lattice, where u = z^^"^ = e~^^ , Uc,sq = 3 — 2^/^, and Uc,t = 1/3- For the honeycomb, 
3 ■ 12^, and 4 ■ 8^ lattices we prove an exact symmetry of the reduced partition functions, 
Zr{z, — yu) = Zr{—z,^). This proves that the zeros of Z for these lattices lie on = 1 for 
— 1 < 2; < as well as the Yang-Lee range < 2; < 1. Finally, we report some new results 
on the patterns of zeros for values of m or 2; outside these ranges. 
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The Ising model has long served as a simple prototype of a statistical mechanical system 
which (for d greater than the lower critical dimensionality ci^.c.d. = 1) undergoes a second- 
order phase transition with associated spontaneous symmetry breaking and long range mag- 
netic order. The general zero-field (spin 1/2) Ising model on a lattice A at temperature T 
and external magnetic field H is defined by the partition function 

Z=T. e-^^ (1) 

{o-„} 

with the Hamiltonian 

= — X! ^nJnn'O'n' — H'^an (2) 

<nn'> n 

where cr^ = ±1 are the spin variables on each site n of the lattice, Jnn' is the spin-spin 
exchange constant, and the units are defined such that the magnetic moment which would 
multiply H an is unity. We shall prove a general result for the case of arbitrary J„„/ 
(connecting any two sites n and n') and then concentrate on the usual nearest-neighbor 
model with uniform coupling Jnn' = JSn,n'±ej, where ej is a lattice vector. We use the 
standard notation P = {ksT)'^, K = (3J, h = jSH, z = e~'^^ , u = = e~^^ , and /i = e"^''. 
The reduced free energy (per site) is / = —(3F = lim^r^oo -^"^ lu 2' in the thermodynamic 
limit, where N denotes the number of sites on the lattice. For the uniform case, we can write 

Z = e«''/2)^+'')^Z, (3) 

where q denotes the coordination number of the lattice. For fixed N, Zr is then a polynomial 
in fi and either u (if q is even) or z (if q is odd). The 2D Ising model has the appealing 
feature of exact solvability; the zero-field free energy / and spontaneous magnetization M 
were first derived by Onsager and Yang, respectively ^ (both for the square lattice; these 
solutions were later generalized to other 2D lattices). 

Yang and Lee pioneered a very interesting line of research in which one studies the model 
with the external magnetic field generalized from real to complex values [0, ^. For an 
arbitrary (finite as well as infinite) lattice with ferromagnetic spin-spin couplings Jnn' > 
(and physical temperature < P < oo), Yang and Lee proved a classic theorem stating that 
the zeros of the partition function lie on the unit circle = 1 in the fi plane and pinch the 
positive real fi axis as the temperature T decreases through the critical temperature Tc 0, Q 
(see also Ref. |0.) In the thermodynamic limit, these zeros become dense and determine the 
continuous locus of points in the fi plane where the free energy, for fixed T, is non-analytic. 
Henceforth, we shall concentrate on the thermodynamic limit, and on results from finite 
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graphs which can give insight into this hmit.Q An elementary property of (|T]) is its invariance 
under the simultaneous transformations cr„ —>■ —cTn, h —>■ —h. This implies that the phase 
diagram in the /i plane is invariant under the inversion 

^ (4) 

with the obvious sign flip M{fi) = — M(l//i). With no loss of generality, one may therefore 
restrict one's considerations just to the interior and rim of the unit circle = 1 in the fi 
plane. 

A complementary complexification is to consider the model in zero field with the tem- 
perature generalized to complex values. In this case, one is interested in the phase diagram 
in the complex u plane (complex z plane for odd q). A number of results in this area, both 



exact and from series analyses, have been obtained for the Ising model with spin 1/2 [3-||T3 
(and with higher spin |Tl|, |l5l). The next logical step is to consider the model with both 
temperature and magnetic field generalized to complex values. Perhaps the earliest step in 
this direction was the exact solution by Lee and Yang 0], of the square-lattice Ising model 
for a particular manifold of pure imaginary values of external field, H = i{7r/2)kBT, i.e. 



h = in /2 {jj, = —1). Recently, we have studied this area further |jT6|, [Tj 

In the present note we report some new results on complex-field (Yang-Lee) zeros of 
the Ising model partition function. We first prove a general result. Consider the Ising 
model with uniform ferromagnetic spin couplings. For d > de,c.d. = 1, the model has a 
low-temperature (i.e., small-w for even q, small-2; for odd q) series expansion with a finite 
radius of convergence. This fact allows one to carry out an analytic continuation from 
u ~ 0^ through m = to negative real values of u in the vicinity of the origin (for odd q, 
a continuation from z ^ 0^ to negative real values of z near z = 0). The properties of the 
model are continuous under this analytic continuation; in particular, the zeros of Z remain 
on the unit circle = 1.0 As we shall now show, the actual interval in negative u or z 
in which the zeros remain on \fi\ = 1 depends on the lattice. We specialize to (i) uniform 

^ Several further questions concerning these zeros have never been answered exactly. These include, for 
the case of lattices of dimensionality d > 2 (i) the density of zeros, and (ii) the location of the complex- 
conjugate endpoints of the distribution for T > Tc- Although we shall not need this here, we mention that 
there has been intensive study of these zeros and generalizations since the original works H, ||; some work 
(on regular lattices without quenched disorder) is listed in Ref. 0. For J„„/ < 0, the zeros do not lie on the 
unit circle |/Lt| = 1 |^. For the ID uniform- J case, an exact solution by Yang and Lee ||] shows that these 
lie on the negative real n axis; for d> 2, there arc no exact results on the location of the zeros for J„„' < 0. 
A numerical study was carried out in Ref. ||] . 

^We report results for real u here. We have also calculated zeros for complex u; in this case, the zeros in 
the n plane are not symmetric under complex conjugation. For example, for the square lattice, we find zeros 
which lie along (i) a spiral curve for u = i, and (ii) two linked, translated spiral-like curves for u = e*'^/^. 
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nearest-neighbor couplings Jnn' = J8n±ei,n' and (ii) the thermodynamic hmit (which, as 
usual, can be probed by using sufficiently large finite lattices). 

We consider the square lattice first. We found previously that for = — 1, the locus 
of points in complex temperature (the u plane) across which / is non-analytic is the union 
of the unit circle and a finite line-segment: 

{U = n {l/We < M < Me} (5) 

where < < 27r and the inner endpoint. Me, of the line segment is given by 



Uf, 



-(3-2^/2) ^ _0,i7i573___ (6) 



Note that Me(/i = —1) = —Uc^sq{,l^ = 1), where Uc,sq is the critical point separating the Z2- 
symmetric, paramagnetic (PM) phase and the broken-symmetry, ferromagnetic (FM) phase 
of the H = model on the square lattice. Now we switch our view back from the non- 
analyticities of / in the u plane as a function of /i to the non-analyticities in the /z plane as 
a function of u. From our previous result (^), we know that as u moves leftward from the 
origin along the negative real u axis, there will be a point of non-analyticity when it reaches 
the value Me, and this point will occur, in the /i plane at = —1. 

To investigate the situation in the interval Ue < u < 0, we use two specific methods: 
(i) exact calculation of the partition function and numerical evaluation of the corresponding 
zeros in /i on finite lattices; and (ii) analysis of low-temperature, high-field series expansions. 
All of our finite-lattice calculations for We < < yield zeros in /i which, to within numerical 
accuracy (~ O(10~^) lie on the unit circle = 1. In Fig. 1(a) we show a plot of zeros of 
Z in fi for an 8 x 8 lattice with u = —0.1. For this and the other figures, periodic bound- 
ary conditions are used. We have also made calculations with helical and open boundary 
conditions, and obtain the same conclusions. 

As u decreases through the value Ue, the situation becomes more comphcated. We find 
that the zeros no longer all lie on the circle |/i| = 1. For u slightly more negative than 
Ue, we find that some zeros lie on the negative real fi axis (in a manner symmetric under 
fi ^ 1/fi, as implied by eq. @). This is illustrated by Fig. 1(b), for u = —0.25. In the 
thermodynamic limit, these presumably merge to form a line segment which originates at 
the point fi = —1 when u = Ue and spreads outward from this point as u moves to larger 
negative values. Further structure may develop around fi = —1 (see Fig. 1(c) for u = —0.5). 
As u becomes more negative, zeros appear elsewhere, e.g. on the arcs evident in Fig. 1(c). 

We have confirmed the locations of these singularities by analyzing low-temperature, 
high-field series expansions. With the definitions 

f=iq/2)K + h + friu,fi) (7) 
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and fr = lim7v->oo N ^ \nZr, these expansions can be expressed as a small-/! series, 

oo 

fr = Y. Ln{u)li^ (8) 
n=l 

or a small-M series, 

oo 

/, = ^V'n(/xK (9) 
n=l 

where L and ip are polynomials in u and /x, respectively. For lattices with odd g, similar 
series hold with u replaced by z. Comparisons of series analyses with exact results, e.g. for 
the 2D Ising model with = (/i = 1) and h = 111/2 {fi = —1; see Ref. show that 

series of reasonable lengths yield very accurate determinations of the positions of singular 
points. The series (H) has been computed to order /i^^ for the square (sq) lattice and 



to order /i for the triangular lattice 0, 19], to be used below. The series @ has been 



calculated to 0(m^^) for the square lattice [^] and to 0(m^^) for the triangular lattice 



For a given value of m, we use the small- /i series to compute the magnetization M and fit 
this to a leading singularity 

M,,„, (10) 

where /i^ denotes a generic singular point (depending on u). The symmetry (^ automatically 
implies that M then has the singularity 

Msing ~ [(1 - /^^s)(l - /is/i)]'/'' (11) 



As in our earlier papers |[T^, |T6|, , we use dlog Fade and (first-order, unbiased) differential 
approximants for our study. Details of our methods are discussed in these papers. We have 
checked our results by analyzing the small-w series (^), to compute M, fitting it to a leading 
singularity of the form 

Ms^ng ~ (1 - u/Usf' (12) 

where Us is a generic singular point (depending on /i). That is, we obtain a series of pairs of 
singular points (/i, u)^; these are singular points in the /i plane for a given u or, equivalently, 
singular points in the u plane for a given /i. The results for the singular point Ms(/i) are 
in excellent agreement with the values of fis{u) obtained from the small-/! series. As we 
decrease u from through real values, we first find a firm indication of a singularity when u 
passes through Ue] this occurs at fi = —1. This singularity moves inward toward the origin 
in the /i plane as u moves to the left of u^- Some typical values are shown in Table I. 

These results are in very good agreement with our findings from the calculation of zeros of 
Z, and suggest that the singular points fis{u) and l/fis{u) are the inner and outer endpoints 
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u 






-1/4 


-0.615(10) 


-0.20(2) 


-1/3 


-0.400(5) 


-0.18(3) 


-1/2 


-0.219(2) 


-0.19(1) 


-2/3 


-0.1407(2) 


-0.20(1) 


-1 


-0.0735(1) 


-0.20(1) 


-3/2 


-0.0372(5) 


-0.20(1) 



Table 1: Values of fig and l/Sg for various u from our analysis of the small-/! series for M 
from (P). For u = u^, the location of the singular point is known exactly as fig = ~1- 



fi 


Us 




-0.15 


-0.6400(5) 


-0.185(15) 


-0.2 


-0.5314(5) 


-0.19(2) 


-0.4 


-0.3345(4) 


-0.19(2) 


-0.6 


-0.2529(3) 


-0.20(3) 


-0.8 


-0.2058(4) 


-0.21(2) 


-1 


= -0.171573... 


-1/8 



Table 2: Values of Ug and (3s for various n from our analysis of the small-w series for M from 
(1). The entries on the last line are exact {ug is given in eq. (|D). 
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of (what in the thermodynamic hmit becomes a dense) hne segment of zeros of Z (i.e., 
non-analyticities of /). 

As Table 1 shows, we find that M has a divergent singularity at /is (and hence, by (0), 
also at 1/ ^is)- This is an analogue in the /z plane of the general result in the u plane which we 
have found in our previous work on complex-T singularities |jl2|, |15|, Q that M diverges, as 
a function of m, at the endpoints of arcs or line segments of singularities protruding into the 
FM phase. In particular, for /i = —1, the exact result for M{u) exhibits divergences at 
u = Ue and u = l/ue, the endpoints of the line segment in (§), with /3e = —1/8, while from 
the exact result 0] for / one can extract the specific heat exponent a'^ = 1. Our series 
analysis strongly suggested the exact value 7^ = 5/4 at this point for the susceptibility 
exponent (so that a'^ + 2/?e + 7e = 2). Although we have found violations of scaling relations 
for complex-T singularities, we note that, provided the usual relations a' + l3{S + 1) = 2 and 
7' = P{6 — 1) hold at M = Ue, it follows that l/6e = —1/9 at (/i, m) = (— 1,'Ue)- The fact 
that this differs from the value 1/5 ~ —0.2 for the (/i,M) entries in Table 1 is not surprising. 



because the point /i = — 1 is quite special, being related by a simple transformation p2|, |T6 



to the zero-field model, i.e. to the point /i = 1. Thus, just as a small nonzero value of the 
H fundamentally changes the properties of the zero field model {H is a relevant parameter), 
i.e., the singular properties change abruptly when /x moves slightly away from /i = 1, so also, 
one expects that a small change in n away from /i = — 1 will have a similarly abrupt effect 
on the singular properties. 

As the results in Table 2 show, M is also divergent at the singular values of Ug corre- 
sponding to the /X values listed there. The corresponding exponent /? ~ —0.2, and, again, 
it is not surprising that this differs from the exactly known value /3e = — 1/8 at m = Ue, 
/i = — 1, for the same reason as given above. 

For large negative u, we can observe that in the polynomial L„ at n'th order in /i in (H), 
the dominant contribution is given by the term of highest power in u, oc u'^^'^. Now consider 
the limit |m| — 00, /i — > with 

X = u^^V (13) 

fixed. Then the double series (^ for fr reduces to a series in the single variable x. We have 
analyzed this series to determine the singular point Xg- For a given large negative value of 
u, we can then extract the asymptotic value of the singularity fis- Clearly, as u — > —00, fis 
approaches the origin like 

where a is a constant (and g = 4 and 6 for the sq and tri lattices). 

We also observe from our calculations of zeros that for small negative u (e.g. Fig. 1), 
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the density g{6) of zeros is consistent with being constant on the circle. Although g{6,u) is 
not known exactly even in the Yang-Lee region < n < 1 or the subinterval of this region 
lying in the FM phase for the square lattice, viz., , < m < Mc (in which the zeros cover the 
entire circle) , we have found by explicit calculation on finite lattices that for u in this latter 
interval that as m — O"*", the density g{6) is again consistent with being constant. 

We have also calculated zeros for real u < 1. These fall in more complicated patterns. 
Here our emphasis is on determining the boundaries of the intervals in negative u where 
the zeros still lie on the unit circle = 1 and establishing the nature of the singularities 
(which we show are line segments on the real /i axis) which first appear when these zeros 
start moving off the circle. 

Next, we carry out an analogous study for the triangular lattice of the size of the interval 
along the negative u axis for which the zeros stay exactly on the circle = 1. We bring 
to bear our knowledge of the complex-temperature phase diagram: for /i = (/i = 1), the 
continuous locus of points where / is non-analytic is the union of the circle and semi-infinite 
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{u = -1/3 + (2/3)e^^} n {-oo <u< -1/3} (15) 

where < < 27r. Recall that the physical critical point is Mc,t = 1/3. For /i = —1, the 
non-analytic points in the u plane consist of a circular arc and semi-infinite line segment ||lj 



{u = (l/2)(-l + e'^) ; 9ce < \0\ < vr} n {-oo <u< -1/2} (16) 

where Uce = e^^"^ = (1/9)(~1 + 2^/^i). From explicit calculations of zeros in /i and series 
analysis, we find that for —1/3 < u < 0, the zeros of Z continue to lie on the circle = 1. 
(The positions are exactly on this circle to within the numerical accuracy of ~ O(10~^).) 
Fig. 2(a) is a plot of zeros for u = —0.1 on a triangular lattice^. As u decreases through —1/3 
moving to larger negative values, some zeros appear on the positive real fi axis, starting at 
/i = 1 and spreading outward along the positive real axis from this point. Fig. 2(b) shows 
the zeros for u = —0.4. It is plausible that in the thermodynamic limit, these merge to form 
line segments. We have used the small-/! and small-w series to determine the position of the 
inner endpoint, fie,rhs, of this line segment (and hence, by the /i — > l//i symmetry, also the 
outer endpoint). We find results in agreement with the zeros calculated on the finite lattice. 
Here, Ihs, rhs denote "left-, right-hand side". As u reaches —1/2, our exact results in Ref. 



T7| show that there is a new set of singularities in the fi plane first appearing at /i = — 1. 



This singularity fie,ihs, moves inward toward the origin along the negative real fi axis (and 

■^The triangular lattice can be represented by a square lattice with an interaction between, e.g., the spins 
on the lower-left and upper-right corners of each square. We use a 6 x 6 lattice of this type with periodic 
boundary conditions. 
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its inverse moves outward) as u decreases past —1/2 toward larger negative values. We find 
that the small-/i series (H), to the order calculated, is not sensitive to fie,rhs] however, the 
small-M series (P) does allow a rough determination of its value. 

Finally, we prove a theorem concerning the reduced free energy for the honeycomb (hexag- 
onal) lattice and for two heteropolygonal lattices: the 3 ■ 12^ and 4 ■ 8^ lattices. (For notation, 
we refer the reader to our previous paper on heteropolygonal lattices [|13|-) These lattices 
all have odd coordination number (g = 3). 

Theorem: For lattices with odd coordination number q, 

Zr{z, -fl) = Zr{-Z, /i) (17) 

This is proved as follows. Now /i — /i corresponds to h h + in /2. Using the identity 

gj(7r/2)<T„ ^ ^^^^ h.me 




Since q is odd, we may replace o"„ = cr^ on each site. We can then associate each factor in o"^ 
with one of the q bonds adjacent to this site. The product of cr^ over sites is thus re-expressed 
as a product of pairs o-„cr„/ over bonds. We next use the identity 6**^'^/^^°'"°'"' = ianCn' to re- 
express this product over bonds in terms of an exponential with shifted coupling K: 

Z(ir,/i + z7r/2) =z(2-'?W2)^gi^E„„,-na„,+/.E„-" =,(2-a)A^/2^(;r^ (19) 

where 

^ = K + Y (20) 

The shift (pO]) takes z = e~^'^ to —z; using this together with (^ then yields the theorem. 

A consequence of this theorem is that for a lattice with odd q, the zeros of the partition 
function in the fi plane lie on the unit circle = 1 not just for the Yang-Lee range < z < 1, 
but for the larger range 

-l<z<l (21) 

In the thermodynamic limit, for Zc < z < 1 or — 1 < z < —Zc, (where Zc = e'"^^" is the 
respective critical point in the z variable for each such lattice) these zeros form a circular arc 
which does not completely enclose the origin in the fi plane, allowing an analytic contination 
from the interior of the unit circle to its exterior. As z decreases through Zc or —z increases 
through —Zc along the real z axis, the arc closes to form a complete circle which prevents an 
analytic continuation from the interior to the exterior of the unit circle = 1. In passing, we 
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recall that for the honeycomb lattice, = 2 — ^/3 (see Ref. [0 for the complex-temperature 
phase diagrams of the Zc values for the 3 ■ 12^ and 4 ■ 8^ lattices). 

This research was supported in part by the NSF grant PHY-93-09888. 

Figure Captions 

Fig. 1. Zeros of Z in the fi plane for the Ising model on a square lattice of size 8x8 for 
u = (a) -0.1, (b) -0.25, (c) -0.5. 

Fig. 2. Zeros of Z in the /i plane for the Ising model on a triangular lattice of size 6x6 
for (a) u = -0.1, (b) u = -0.4. 
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